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We study the edge states in a two dimensional electron gas with a transverse magnetic field and Rashba 
spin-orbit coupling. In the bulk, the interplay between the external field perpendicular to the gas plane and the 
spin-orbit coupling leads to two branches of states that, within the same energy window, have different cyclotron 
radii. For the edge states, surface reflection generates hybrid states with the two cyclotron radii. We analyze the 
spectrum and spin structure of these states and present a semiclassical picture of them. 
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I. INTRODUCTION 

Since the seminal spin-transistor proposal by Datta and 
Das^ it has been recognized that the spin-orbit interaction 
may be a useful tool to manipulate and control the spin degree 
of freedom of the charge carriers. This opens novel opportu- 
nities for the developing field of spintronics.^ The challenging 
task of building spin devices based purely on semiconducting 
technology requires to inject, control and detect spin polarized 
currents without using strong magnetic fields. For this pur- 
pose, the spin-orbit coupling may be a useful intrinsic effect 
that links currents, spins and external fields. During the last 
years a number of theoretical and experimental papers were 
devoted to study the effect of spin-orbit coupling on the elec- 
tronic and magnetotransport properties of two dimensional 
electron gases (2DEG).^^^^^^^^W^^^W^ This is endorsed by 
the fact that, in some of the semiconducting heterostructures 
used to confine the electron or hole gas, the spin orbit interac- 
tion is large. Moreover it may be varied by changing carrier 
densities or gating with external electric fields. ^"^'^^ 

In many transport experiments in 2DEG with a transverse 
magnetic field, including quantum Hall effecU^ and transverse 
magnetic focusing, ^^'^^'^^ edge states play a central role. To 
our knowledge, a detailed analysis of the effect of the spin- 
orbit coupling on the edge states has not been done yet. In 
this paper we present a theory for edge states in 2DEG with 
transverse magnetic fields and a Rashba term describing the 
spin orbit interaction.^^ 

First we focus on the quantum mechanical solution. By 
numerical diagonalization of the Hamiltonian in a truncated 
Hilbert space we calculate the energy spectrum and the wave- 
functions that, as we show below, present an intricate struc- 
ture. Then we resort to a semiclassical analysis to interpret 
and illustrate the nature of the edge states in the high energy 
or low field limit. 

Our starting point is a 2DEG with Rashba coupling and an 
external magnetic field B perpendicular to the plane contain- 
ing the electron gas: 

where m* is the effective mass of the carriers, P^^ = -\- 
{e/c)Arj, with Pn and Arj being the r^-component of the mo- 
mentum and vector potential respectively, a is the Rashba 



coupling parameter, are the Pauli matrices and g is the gy- 
romagnetic factor. The last term V{x) is the lateral confin- 
ing potential. For simplicity, from hereon we consider a hard 
wall potential that confines the electrons in the transverse x- 
direction: V{x) = for < x < L and infinite otherwise. 



II. THE QUANTUM SOLUTION 

In the geometry where electrons are confined in the x- 
direction, it is convenient to use the Landau gauge A = 
(0, xB^ 0) and write the wavefunction in the form: 



^{x,y) = e'^y(p{x) 



(2) 



with the function (p{x) expanded in the basis set of the infinite 
potential well 



(p{x)-- 



The Schrodinger equation = leads to the follow- 
ing equations for the spinors 




Y}^Pl^-Glm)<J--{Fi^ + Glra)(J+-MiJ^ 
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with Mim, and Gim proportional to the matrix elements 
of the operators (x — xq)^, (x — xq) and d/dx respectively. 
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Here, c^c = e |5| /m*c is the cyclotron frequency, — (gx ± 
icFx)/2 and xq = —Hkc/eB. We solve these equations in a 
truncated Hilbert space disregarding the highest energy states. 
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FIG. 1: Probability and spin densities for the 1^* and 2"^"^ levels (left 
panel) and 3^"^ and 4*^ levels (right panel) respectively for xo = 0, 
m* =0.068me, a^lOmeVnm, L = 600nm and B = 2.5T. 




40 80 120 



X [nm] 

FIG. 2: Charge and spin density current for three values of the Fermi 
energy (left panels). The Sz (thick line) and (dotted line) current 
densities are measured in units of h/2e. Charge current density is 
indicated with a thin line. Note than in (b) and (c) the current density 
is multiplied by a numerical factor indicated in the figure. In (d) 
the energy of the first levels versus xo are shown. The three values 
of the Fermi energy used in (a), (b) and (c) are indicated. In (e) 
the semiclassical orbit is shown with the velocity and spin direction 
indicated by arrows. Parameters as in figure 1. 



Typically we take a matrix Hamiltonian of dimension of a few 
hundreds and keep the first thirty states. In all cases the width 
of the sample L is taken large enough to have the cyclotron 
radius Tc smaller than L/2. The right and left edge states are 
then well separated in real space. For xq L/2 the states are 
equal to the bulk states, except for exponential corrections. 
The wave functions and the energy spectrum reproduces the 
known results: in the bulk the spin-orbit coupling mixes the 
two spin components and there are two branches of states with 
energies given b}iS2iSi 



E^ = hu;cnT^Ei^^-j 2n (6) 

with n > 1 and a single state (n = 0) with energy Eq = 
huJc/2 — giiBB/2. The corresponding eigenfunctions for n > 
1 are 

and 

Here Ly is the length of the sample in the ^/-direction, 0n(^) 
is the harmonic oscillator wavefunction centered at the coor- 
dinate Xo, An = l -\- D'^ and 




The wavefunction of the state with n = is 



^0,/c(^,y) = 




In the bulk (xq ~ L/2), the ground state has the spin 
along the z-direction. In the excited states the spin is tilted 
with an expectation value of its z-component (a^) = ±(1 — 
\Dnf)/An that decreases as a and n increase. The condi- 
tion Eq « 2{a/lcf' , that is equivalent to (a^) for all 
n > 0, is referred as the weak field condition iSi For large 
enough n, when Eq << 2n (a//c)^, the spin-orbit dominates 
and {az) ~ 0. 

For high fields or low electron density, the physical proper- 
ties of the system are dominated by the states with low quan- 
tum number n. We first consider this case and present the 
results for the first few Landau levels. 

As the momentum k parallel to the edge varies, the cen- 
ter of gravity xq of the wavefunctions changes and as it ap- 
proaches the sample edge,the effect of the confining poten- 
tial becomes important generating the /c-dependent dispersion 
of the energy levels. ^^The interplay of the spin-orbit coupling 
and the confining potential produce a tilting of the spin for all 
the edge states. Edge states probability densities |^|^ and the 
corresponding spin densities Si = (fi/2)^Vi^ are shown in 
Fig [l] for xo = . The spin is predominantly in the xy-plane 
even for the lowest energy edge states and the sign of the spin 
densities alternates as the energy increases. The current car- 
ried by these states is then polarized and for the parameters of 
the figure the polarization is determined by the Rashba cou- 
pling. We use the charge and spin currents operators defined 
as: Jy = —ey for the charge current and J^^ = hyax/2 and 



Jy^ = h{yGz -^cFzij)/^ for the spin currents. ^ In these expres 
sions the velocity operator in the ^-direction is given by 
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(11) 



The total current densities, defined as jy{x) = 

^occ^n'^^^->y)'^y^n{^->y) where the sum runs over 
all occupied states, are shown in Fig|2| for different values 
of the Fermi energy. The charge and jy^ current densities 
are confined at the sample edge indicating that they are due 
to the edge states. Conversely the current density has a 
non-zero value inside the sample. The origin of this current 
can be understood in terms of the simple semiclassical picture 
shown in Fig|2t: electrons moving in the positive (negative) 
^/-direction have a positive (negative) projection of the spin 
along the x-axis. Since the spin is not conserved, these 
currents do not necessarily produce spin accumulation in 
samples with constrictions or edges perpendicular to the 
current direction. 

Let us now consider the low field case where many Landau 
levels are below the Fermi energy. The energy spectrum as 
a function of xq is presented in Fig|3l For the bulk states, 
the typical energy splitting of the two branches (+) and (— ) 
is different (see Eq.©), leading to a beat in the total energy 
spectrum. Within the same energy interval, the two branches 
have different quantum number n and consequently different 
cyclotron radius Vc- We take^^ 



--2{^^\{x-x,f\^i) 



(12) 



that for large n gives 2n{h/m'^u;c). According to equa- 
tion (|6|), in this limit states with approximately the same en- 
ergy belonging to different branches have cyclotron radius dif- 
fering in Arc — 2a/huJc- Note that the radius difference in 
this large n approximation does not depend on n. 

For xq < Tc the effect of the confining potential becomes 
important and the two bulk branches mix leading to edge 
states that combine the two cyclotron radii. This mixing is 
apparent from the energy spectrum that presents level anti- 
crossings as shown in the inset of Fig|3l 

The behavior of the levels 23 and 24 is illustrated in Fig|^ 
In the top panel the figure their probability densities for xq = 
I// 2 are shown. These states correspond to a state of the (+) 
branch with n = 13 and a state of the (— ) branch with n = 10 
respectively. The (+) branch state radius is larger than the (— ) 
branch one as can be inferred from the figure. We can follow 
the evolution of these states as xq changes from xq = L/2 
to a negative value. A contour plot illustrating this evolution 
in shown in the central panel of Fig|4| For xq ~ Tc a sud- 
den change in the wave function spatial extension is observed. 
States belonging to the bulk (+) branch shrink for xq ~ Tc 
due to the mixing with the (— ) branch states. Conversely, 
states belonging to the bulk (— ) branch expand for xq ^ Vc- 
This sudden change in the wave function extension is of the 
order of Arc. The lower panel of Fig|4| shows the probability 
densities of the two levels for a negative value of xq. 
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FIG. 3: Energy spectrum versus center coordinate xq for a system 
with m* 0.068me, a lOmeVnm, L GOOnm and B 0.5T. 
Inset: detail of the anticrossing of levels 23 and 24. 
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FIG. 4: Evolution of the wavefunctions of levels 23 and 24 with the 
center coordinate xq. Parameters as in Figure 1. Upper panel: prob- 
ability densities of the levels 23 and 24 for xo = 300nm. Central 
panel: density plot of the probability densities of the two levels ver- 
sus Xq. Lower panel: probability densities for xq = — lOOnm. 



For Of = the number of nodes of the wavefunction of a 
given level is conserved as xq changes. With spin-orbit cou- 
pling, the anticrossing of energy levels is an indication that 
the wavefunctions change in character as xq changes. Then, 
for a given energy level the number of maxima of the proba- 
bility density is no longer conserved as it is shown in Fig|4| 
It is also interesting to analyze the spin densities associated to 
these states. The spin structure of the edge states for xq =0 
is shown in Fig|5l The spin densities Sx and Sz show an intri- 
cate behavior due to the beating of two contributions. This is 
a consequence of the mixing of states that are (+) and (— ) in 



4 



character. As we discuss below the semiclassical analysis also 
unveils that edge states are formed by combining states with 
different radii and different spin projections. 



III. THE SEMICLASSICAL SOLUTION 

A. Bulk States 

In a recent approach,^"^ both the orbital and the spin degrees 
of freedom have been treated semiclassicaly in an extended 
phase space. The spin coherent state is defined as 



\z)- 



i) 



(13) 



where z is a c-number. A unit vector associated with the clas- 
sical spin is defined in terms of the coherent state \z), which 
for spin | reads 



n={z\a\z) , 
with components determined by z 

2z* 



ni + in2-- 



(14) 



(15) 



and nf-^ n^-^ ^3 = 1- 

The classical phase-space symbol H( q, p, n) of the Hamil- 
tonian defined in Eq.Jl]) iJ 



24 



H(q,p,n)=Ho(q,p)+ -n-C(q,p) (16) 

where the first term is the classical Hamiltonian without spin 
orbit coupling and 



C(q, p) = {^{Py 
The equations of motion are: 



e 2a 

-Bx),-—p,,0). (17) 



dH 
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dH 
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n = C X n 



(18) 



whose solutions represent the classical orbits in the ex- 
tended phase space. In bulk there is an additional constant 
of motion besides the energy and the system is classically 
integrable. We are interested in the periodic solutions of 
Ea.fTst from which the action integral can be computed in 
order to apply an EBK quantization scheme. We propose 
q = r(cos cjt, sin cot) and replace it in Ea.fTsb. For t = the 
initial conditions Px{0)=Py{0) = SLve impossed leading to 



eB 
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FIG. 5: Probability densities and the corresponding spin densities for 
levels 23 and 24 and xo = 0. Parameters as in Fig|3l 



where we have defined A = u; — cOc- The normalization condi- 
tion for the spin components gives an algebraic equation of or- 
der fourth in A that to leading order in h (semiclassical limit) 
gives 



hr 



(20) 



Given the two frequencies co = cOc^c^/hr one must go one 
step further to find an explicit expression for the cyclotron 
radius. Replacing Eq.(IT91 evaluated a t = in the equation for 
the energy conservation H = E, and taking into account the 
value of A obtained in (l20ll . we find 




(21) 



Therefore for a given energy E the periodic solutions re- 
sult in two orbits of radii r±, frequencies uj± = ujc T 
and opposite values of the spin respectively. The cyclotron 
radii difference is Ar = r+ — r_ = 2a/huJc in exact cor- 
respondence with the quantum mechanical estimate that we 
obtained for large n. In these two orbits, with different radii 
and frequencies, the electron has the same velocity; that is 



r+c.+ =r_c._ = v'(t) 

Once we know the periodic solutions in extended phase 
space we need to compute the action integral /. Following 
Ref. 12411 . the action can be expressed as: 



-L 



pq + -(^1^2 - ^2^i) ) dt (22) 



eTT ^ 9 h 



with T± = 2t: /uo±. The EBK quantization rule sets 
/ 
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- 7 = 27m 



(23) 
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FIG. 6: Skipping orbit used in the semiclassical calculations drawn 
up to the first specular reflection (first period of the x motion). The 
values of x+ and x- have been enlarged in order to visualize the 
angle 77 (see text for details). 



where n is an integer and 7 is the sum of the phase shifts 
acquired at the turning points of the motion^^. For the bulk so- 
lutions the turning points are two caustics giving a total phase 
shift 7 = —TT. Replacing Eq.iEJ in Eq.(l22ll we finally obtain 



E =nhu;c T -j-V2n , 



(24) 



which is the quantum spectrum for the bulk states (neglecting 
the zero point energy)^^. With the notation we emphasize that 
for a given quantum index n, the (energy associated to 
r+) is lower than E~ (energy associated to r_). 



B. Edge States 

As we mentioned above, for xo ~ Tc the bulk (+) and (— ) 
branches mix, and the wavefunctions spatial extension present 
a remarkable change. The avoided level crossings structure 
observed around xq ~ Tc in the energy spectrum is the finger- 
print of this behavior. The semiclassical image that we pro- 
pose consists of a skipping orbit formed by a series of trans- 
lated circular arcs of radii r+ and r_ and centers x+ and X-in 
the X direction respectively (see Fig|6ll. The center coordinate 
y of these circular arcs changes at each specular reflection. In 
Fig |6| this primitive orbit is plotted for a complete period of 
the X motion. The fact that the reflection at the boundary is 
specular is guaranteed by the conservation of the modulus of 
the velocity ( r+cc;+ = r-UJ-) and can be cast in the form 



cos T] = = = C 



(25) 



being r] the angle depicted in Fig|6| 

As it can be inferred from the semiclassical solutions ob- 
tained in Eq.iO, to lowest order in h the in plane spin com- 
ponents (ni and 712) of the orbit r+ have opposite signs than 
those of the orbit r_ and in both cases is 713 = 0. Therefore 
the spin conservation is guaranteed at each specular reflection 
of the skipping orbit with the boundary if x+ and X- ~ 0. 
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FIG. 7: Edge states energy spectrum obtained from the semiclassical 
approach Ea.dZTl (thick solid lines) together with the exact quantum 
results. Around xo = both results are almost indistinguishable even 
in the low energy region. Inset: detail of the comparisson for xo = 
and intermediate energies. 



Our goal is to perform a semiclassical quantization employ- 
ing this classical skipping orbit. The semiclassical approach 
is fully justified for angles r] it/2 and we will obtain the 
energy spectrum and the dispersion relation quite accurately 
around xq ~ 0. 

We proceed analogously to the previous section, taking now 
into account that the orbital motion projected on the x axis is 
periodic with a period Tg = Ti + T2 = 2(7r — 7^)/cc;+ + 2(7r — 
r])/uo-. For the sake of clarity we divide the action integral 
in two terms I = Io + • The action associated to the orbital 
motion is 
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For the spin degrees of freedom the action integral is straight- 
forward to evaluate and gives Is = 2h{iT — r]) . For the skipping 
orbit the phase shift is 7 = tt, due to the fact that we have to 
consider now for each period of motion two bounces with the 
boundary and also two caustics. Replacing the obtained val- 
ues for I = Io -\- Is and 7 into Eq.(l22ll we finally obtain 
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(27) 



as a function of the parameter (. 

For x+=x_=0, isC = and the energy levels are 

E{0)=nhiJc-m* (^^y . 



(28) 
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FIG. 8: The numerically integrated semiclassical results. Semiclas- 
sical skipping orbit with the in-plane spin projection indicated by ar- 
rows (upper panel). The three components of n as a function of time 
(lower panels). The energy of the orbit is about 37hiOc- Parameters 
as in Fig. 3. 

To obtain the dispersion relation, we proceed numerically due 
to the fact that the variable ( depends on the energy through 
the cyclotron radii. In order to compare with the quantum me- 
chanical solution for C 7^ one needs to rewrite Eq.dZTl as 
a function of the center of the classical skipping orbit which 
plays the role of xq for the edge states. For ^ ~ we have 
checked that = x+ , = X- or xq = (x+ + x_ ) /2 leads to 
almost the same dispersion relation. In Fig (8] we plot Eq.iIZTI 
after chossing xq = (x+ + x_)/2. Notice that the semiclas- 
sical solution follows quite satisfactory the quantum results 
even in the low energy region and it is almost indistinguish- 
able from them for xq ^ 0. 

We end this section by presenting results of the numeri- 
cal integration of the semiclassical equations Eq.dTsl using a 
fourth order Runge-Kutta algorithm. The results agree with 
the analytical solution obtained to lowest order in h, and they 
are summarized in FiglS] where the skipping orbit with two 
radii is clearly observed. The components ni and n2 of the 
classical vector n show the expected behavior with a contin- 
uous evolution at the bouncing point. The out of plane com- 
ponent ns is small and presents fast changes at the bouncing 



points. This component decreases as the energy increases in 
agreement with the semiclassical assumption that, to lowest 
order in h, predicts ns = 0. 



IV. SUMMARY AND CONCLUSIONS 

We have analyzed the eigenstates and the energy spectrum 
of a 2DEG with spin-orbit coupling in the presence of a per- 
pendicular magnetic field. We focused on the edge states 
that appear when the 2DEG is confined in the transverse x- 
direction by a square well potential. We first discussed the 
low energy states in the high field limit. The rest of our work 
was devoted to study the high energy states (high quantum 
number n). In this regime, the spin-orbit coupling has an im- 
portant effect on the edge states: while for the a = the edge 
states with k = (that corresponds to xq = 0) have an energy 
separation AE = 2hjOc}^ for a / the energy separation is 
hwc^ see Eq. (l28b . As pointed out in section II, the effect of 
the spin-orbit coupling increases with a and with the quantum 
number n, and there is always a high energy regime where the 
spin-orbit coupling dominates. In this high n regime, the en- 
ergy spectrum of the edge states follows Eq. (l28t . 

In the bulk, states with large and small radii are quasi- 
degenerated. The bouncing at the surface mixes them lead- 
ing to hybrid states that combine large and small radii. The 
mixing is evident in the quantum solution where the energy 
spectrum versus xq shows avoided level crossings, a fingertip 
of level mixing. 

The spin texture of these states is also discussed in terms 
of the classical solution. To lowest order in h the spin lies 
in the plane of the 2DEG and its direction is perpendicular to 
the velocity. The relative orientation of the spin with respect 
to the velocity is different along the segments with large and 
small cyclotron radius. 

The picture obtained with the semiclassical approximation 
accounts for the quantum mechanical prediction of a splitted 
transverse focusing peaks^^ that was recently experimentally 
observed in hole gas in GaAs.^^ 
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